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1.  Introduction.  In  this  paper  we  describe  a  systematic  approach  to  the 


dimensional  reduction  of  nonlinear,  monotone  boundary  value  problems.  Our 
work  is  motivated  by  problems  arising  from  considering  "thin"  mechanical 
structures  such  as  beams,  arches,  plates  and  shells,  cf.  [2]  and  [26].  Let 
»  ]0 ;  1[x]  -  e;  e[  be  a  domain  in  R^.  We  present  an  efficient  procedure 
to  produce  approximate  solutions  to  the  above  boundary  value  problems  defined 
over  £2e  .  Led  by  intuition,  we  expect  that  the  solution  can  be  approximated 
well  by  linear  combinations 


j=0 


(x2/e) 


x1  [0 ;  1 ] ,  x2  [-e,e] . 


We  select  the  basis-  or  "Ansatz"  functions  {<|>  }"  q  so  as  to  obtain 

maximal  order  of  convergence  as  e  tends  to  zero,  _  are  found  as 

J  j  =0 

solutions  to  Galerkin  systems  of  ordinary  differential  equations.  We  consider 

three  asymptotic  ranges  of  loads  (infinite,  finite,  and  zero  limit)  and  the 

selection  of  "Ansatz"  functions  is  shown  to  be  independent  thereof.  This  is 

of  course  of  major  importance  for  practical  computations.  Estimates  for  the 

order  of  convergence  as  N  tends  to  infinity  are  established  (as  they  must 

be,  since  e  is  never  infinitesimally  small  in  practical  applications  and 

since  a  required  accuracy  is  often  not  obtainable  by  taking  more  terms  in  the 

r  i  N 

asymptotic  expansion).  The  span  of  I’J'jfj.o  is  shown  t0  be  °Ptimal  in  the 
sense  of  N-widths  in  approximation  theory. 

This  method  differs  from  another  method  of  dimensional  reduction  widely 
used  in  structural  mechanics,  that  of  asymptotic  expansion  in  e  ,  see 
Friedrichs  [8],  Friedrichs  and  Dressier  [9],  Reiss  and  Locke  [28], 

Gol 'denweizer  [11],  and  Ciarlet  and  Destuynder  [6].  The  latter  method  is 

2 


limited  in  its  practical  scope  since  it  deals  with  rough  input  data  by 


introducing  boundary  -  and/or  interior  layer  expansions  which  often  are  diffi¬ 
cult  to  obtain  without  extensive  analysis  dependent  upon  the  particular  case 
at  hand.  This  limits  the  generality  of  the  latter  approach.  We  propose  to 
deal  with  such  nonsmoothness  by  increasing  N  near  the  layer. 

The  aforementioned  characteristics  of  our  approach  are  practically 
feasible  as  witnessed  through  our  computational  experience  in  the  last  section 
of  the  paper. 

We  remark  that  this  work  was  inspired  by  the  papers  of  Vogelius  and 
Babuska  [30,  31,  32]  on  dimensional  reduction  for  linear,  elliptic  boundary 
value  problems. 
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2.  Notation  and  model  problem. 

Let  A  be  an  operator  of  the  form 

(1)  Au(v)  »  J  F( | Du| 2)Du*Dv  dx, 

fle 

where  x  «  (x^Xg)  (  fle  -  ]0,1[x]-e,e[  .  Denote  by  »  [0,1]  x  {+e} 

and  [0,1]  x  { -e }  respectively,  and  by  rE  =  ({o}  x]-e,e[)t»  ( { 1  }*] -e , e[ )  , 
so  that  3fle  =  rj  U  r£_  UTq  disjointly .  Let  Du  =  Ou/3x1 ,  3u/3x2>  and 

|  Du  |  «  (  f  (3u/3x.)2)  1/2  . 

i-1 

Let  for  some  n  (  Z+, 

(2)  F(t)  -  1  +  tn  ,  Vt  £  R+  {0 }  , 

for  the  sake  of  explicitness.  What  is  to  follow  could  be  extended  to  other 
choices  of  F  without  much  imaginative  effort,  cf.  [27]. 

We  consider  the  minimization  problem:  Find  uE  such  that 


u^W1'2^2  (ne)  and  E  (ue)  -  inf  1f2n+2  E(v)  , 

v  *  w  (0) 


where 


..1 ,2n+2  ,  e .  ,,1  ,2n+2,„e..  ~  r  .  -i  , 

W  ^  (fi  )  =  W  (0  )  H  {v  :  v|  £  -  0}  ,  and 

ro 

E(v)  -  B(u)  -  G(v), 

where  the  Frechet  derivative  of  B  at  u  is  Au,  and  where  G  k  (W1!2^*2) 
-  the  dual  space  of  W1 |2l^+2  .  is  defined  for  all  y  £  R  and  all 
8  k  L2n+2(0, 1 )  through 

G(v)  -  Jq  8(x1 )  e1  ^[v(x^ ,e)  +v(x1>-e)]  dx1 


Here,  y  is  introduced  in  order  to  investigate  three  asymptotic  ranges  of 


"loads"  on  rE  .  According  to  whether  p  >  0,  p  =  0,  or  p  <  0  ,we  have  an 

infinite,  finite,  or  zero  "load"  (e1  can  be  viewed  as  a  scaling  constant;  it 

X2 

is  the  Jacobian  of  the  transformation  (x1  ,x^)  -*■  (c,n)-(x1(— )  ). 

The  minimization  problem  (3)  is  formally  equivalent  to  the  homogeneous 
Neumann/Dirichlet  boundary  value  problem 

i 

-D*(F( |Duej2))  Due  =0  in  fle  , 

(4)  u£  =  0  on  Tq  , 

F(|Due|2)  3ue/3v  =  8e1_u  on  rE  , 

where  v  is  the  outward  normal,  and  furthermore,  (3)  is  equivalent  to  finding’ 
ue  i  W1 [o"+2  satisfying 

(5)  v  v  *  w1(q)+2  (fl£)  *  Aue(v)  -  G(v) . 

Because  E  is  strictly  convex  and  lim.  .  E(v)  =  ®,  (3)  and  (5)  have  a 

|v|++® 

unique  solution  uE  w1(o^+2  »  cf.  Visik  [29],  Morrey  [25],  Glowinski  & 

Marocco  [10],  or  Ciarlet  [7];  we  define  |v|  =  (  /  /  |Dv|2n+2dx)  1/2n+2 

fle 

(using  Friedrich's  inequality). 

Finally,  we  scale  the  problem  as  follows 

(6)  ue(x1tx2)  -  ue(x1,x2/e)  ,  (x1,x2)  £  Ge  . 

Denote  by  (5th)  =  (x^x^e)  ,  D£  =  (3/3£,  e  ^/Sn)  ,  and 

0  -  ]0,1[x]  -  1 , 1  [  .  (5)  becomes:  Vv  (  W1^*2  ^ 


(7. a) 


f  F( |D  u£|2)D  u£  •  D  v  e  d£dn 
;  fl  1  e  1  e  e 

=  e  JJ  B(5)[v(5,1)  +  v(C.-1)]  e"yd*i  . 

—  £ 

In  the  remainder  of  the  paper,  we  shall  denote  u  by  u  ,  and  the  left  hand 
side  of  (7)  by  Au(v)  and  the  right  hand  side  by  G(v)  for  any 
u,  v  fc  W1  |2l^+2( fl)  .  Thus  (7)  can  be  written  (find  u  i  W1  (o^+2  ^  such 
that ) 

( 7 -b )  Vv  (  W1(2"+2  (Q)  ,  Au(v)  =  G(v) . 

Also,  w1(q^+2  is  endowed  with  the  norm  |v|  =  (Jn|Dev|2n+2ed^dri)1/2n+2 

(again  with  the  use  of  Friedrich’s  inequality;  here  jft|  =  2  independent  of 

e  ).  Note  that  |v£|  =  |v| 

fl£ 

The  scalar  problem  (4)  is  a  mathematical  model  of  the  physical  models  of 
antiplane  shear  in  finite  elasticity,  see  Gurtin  &  Temam  [12],  Abeyaratne  [1], 
and  Knowles  [16,  17],  and  the  torsion  problem,  see  Kachanov  [15]  and  Langen- 
bach  [20].  We  are  trying  to  express  a  nonlinear  constitutive  (stress-strain) 
law,  so  this  covers  finite  elasticity  only  insofar  as  there  exists  some  in¬ 
herent  scalar  quantity  in  our  problem.  In  [1]  e.g.,  the  stored  energy  func¬ 
tion  W  is  assumed  to  depend  only  on  the  first  invariant  of  the  left  Cauchy- 
Green  strain  tensor:  3  +  | Du |  ,  where  the  inherent  scalar  is  u(x1tx2)  the 
x-j-displacement  of  a  prismical  body  along  the  x^-axis,  see  [1]  sections  [2] 
and  [5].  In  the  absence  of  body  forces,  the  problem  is  described  by 
D(W'(3+|Du|2))Du  =  0  in  n  with  boundary  conditions. 

The  most  direct  interpretation  of  (4)  is  that  of  the  (stationary)  heat 
equation  with  a  nonlinear  thermal  conductivity  (dependent  upon  the  temperature 
gradient) . 
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3.  The  Galerkin  approach  to  dimensional  reduction. 

In  this  section  we  present  some  known  results.  We  will  consider  the 

Galerkin  approach  to  solve  (7)  and  (2)  on  a  proper  subspace  VN  which  is 

r  i  N 

characterized  by  the  Ansatz  functions  ill.}.  _  : 

J  J  =0 


VN  -  {  l  c  (O*  (n):{c  }  S  w1 'ln+2  (0,1)} 

j-0  J  J  J 


More  precisely,  define  uN  to  be  the  solution  in  VM  of 


Aun(v)  =  G(v)  ,  V  v  €  VN, 


where  A  and  G  were  defined  in  the  last  section’s  closing  paragraphs.  We 
note  that  the  problem  (9)  is  equivalent  to  the  minimization  problem  (3) 
restricted  to  VN,  which  had  a  unique  solution  (see  references  mentioned  in 
connection  with  (3),  same  argument). 

The  operator  AsW1^"*2^)  -  (W1  | q"+2(  Jl) )  *  is 

i)  strongly  monotone,  i.e.,  j  x:R+  ^  { 0 }  -*■  R  ,  strictly  increasing 

with  x(°)~0  and  limt_><jox(t)=®,  such  that  (Au-Av)(u-v)  2  x(  |u-v  | )  ||u-v  |  , 

y  u,v  t  W  ^  (fi)  ,  and 

ii)  Lipschitz-continuous  for  bounded  arguments,  i.e.,  V  r>0  J  r(r) 
such  that  if  |u|,|v|  £  r  ,  then  |Au-Av|  £  T(r)  Ju-v|  ;  cf.  Ciarlet  [7], 
sec.  5.3. 

This  is  proved  easily  along  the  lines  there  which  are  based  upon  the  work 
of  Glowinski  &  Marocco  [10]  and  could  be  generalized  to  hold  for  polynomials 
F  that  satisfy  (sF(s2)  -  tF( t2) ) /(s-t)  >  0,  V  s,t>0,  s  *  t,  see  [1 U].  It  is 


also  known  (cf.  Ciarlet  [7]  or  Glowinski  &  Marocco  [10]  that  when  A  satisfies 


these  two  properties,  then  there  exists  C  independent  of  Vjj  such  that 


x(|u-uN|  S  C  inf  |u-vN|  . 


VNtVN 


In  our  case,  3  a>0 :  x:  t  -»at2n+1  ,  Vtfc  0  .  This  is  a  nonlinear  generaliza¬ 
tion  of  Cea's  lemma. 

Since  C  in  the  last  inequality  depends  on  e  ,  we  give  a  proof  of  the 
following . 

Theorem  1  .  For  n  i  N,  let  A:!**' 1 2!J+2( ft)  •*  (W1  ’^2^+2  (£1))  be  defined  as 
above.  Then  there  exist  positive  constants  a  ,  M  ,  C  independent  of  t 
such  that: 


(i) 

(ii) 

(iii) 


(Au-Av)(u-v)  £  a  |u-v| 


2n+2 


| Au-Av {M<  ||v | )2n  +  2en/(n+1  ^ju-vj 


1-u 


1-M 


II  <r  2n  +  1  |  ■  s  r  2n+1 

|u|  SCe  ,  |uN|  <  C  e 


Proof :  (i)  In  Ciarlet  [7],  e.g.,  it  is  shown  that  there  exists  a 


(lx|2n-  |y |2ny)«(x-y)  ^  a  V  x,y  i  1  2,  x  *  y 
i  2n+2 

x-y 


which  yields  (i)  immediately,  (ii)  Also,  in  Ciarlet  [7]  is  given  that 

i 

there  exists  M  : 


S  m'  ,  V  x,y  k  R2  ,  x*y  , 

I  x-y | ( M  +|y| )2n 


which  yields  that 


m 

VA 


s'.  1*1 


$;!*i 

Isji’l 

v,W 


:i*S 

#; 

SW 

>•:*! 


t£& 

>•{ 

>w, 


•wa 

*«$ 

»*;S 


2# 


J Au-Av |  =  sup  |(Au-Av)(w)| 

M*1 


S  sup  {  /  M'|D  u-D  v ( ( D  u|+|D  v | ) 2n { D  wjed£  dn 
|w|s|  a 


+  j  |D£u-Dev|  |D£w|edS  dn} 


S  sup  {m  |u-v|  | w |  (|u|  +  |v|)2n  +  | u-v |  | w |  2e  n  +  1  }  , 
|w|s| 


using  Holder's  inequality  twice  (p^  =  2n+2,  p1  =  2n+2/2n+1 ;  p2  =  2n+1 ,  p2  = 
2n+1/2n)  and  the  fact  (x+y)2n  S  22n~1  (x2n+y2n)  for  x,  y  >  0,  see  Korovkin 
[18]  p.  25.  Thus  (ii)  holds.  For  (iii), 

|u|2n+2  £  Au(u)  =  G(u)  £  |Gj  |  |  , 

with  u  replaceable  with  uM. 


| G |  =  sup  |  /  e1  ug(w( • , 1 )  +  w(  • , -1 ) )d£ | 

M*l  0 


SCe 


since  the  trace  operator  maps  W^2^+2  (P)  continuously  into  L2n+2  (D  with  a 


norm  containing  a  factor  e 


;  (iii)  follows. 


Combining  inequalities  (i),  (ii),  and  (iii)  in  Thm.  1  we  finally  arrive  at  the 


following  pair  of  inequalities  in  (10) 


2nd  -n- 


|u-UN|  S  C  £ 


(2n  +  1  )‘ 


inf  |u-v  |2n+1  ,  for  u  >  0  , 

^  r  *• 


(10) 


v  V 
N  N 


,  ...  (2n+1  )(n+1  )  .  .  -  ,2n+1 

|u-uN|  S  C  e  mf  |u-v  | 

v  f  V 
N  *  N 


for  u  £  0 


where  C  does  not  depend  on  t  . 

Proof  of  (10):  To  arrive  at  the  e-dependence  we  present  a  standard  argument: 


“lu-uN|2n+2  5  (Au-Aun)(u-un) 


(Au-Aun)(u-vn) 


s  |Au-Aun|  |u-vn| 


5  1m(|u|  ♦  |\|)2n  +  2e  n+1  }|u-uN|  |u-vN|  ; 

1-U 
2n + 1 

with  | u |  ,  |uj  £  C  e  ,  the  result  follows. 

Due  to  this  reduction  in  power  in  Cea's  lemma  (by  2n+1 )  it  becomes  of 

1  ? 

interest  to  obtain  a  quasioptimal  bound  in  W  ’  under  further  assumptions. 
This  was  done  by  Rosier  [19].  We  present  the  argument  for  our  case.  The 
final  estimate  is  (13)* 

For  v  (  w! (Q),  let  J v |  =  (  /  | D  v| 2ed£dn) 1 /2 


n+1 


Note,  that 


If  we  now  invoke  the  Mean  Value  Theorem  we  get, 


w  w  £  r\  ’  °°  v  £  w1 ,2n+2 

1 ,w2  (0)  1  ’  (0) 


|  ( Aw1  -Aw2)  (v)  |  =  |  /  [F(|D£W1|2)  D£w1 -F(  |  D£w2|  2)  D^^ve  d£dn| 

(12)  =  |  /  D(F(|6|2)  (D£w1-D£w2)  D£vedC  dn| 

*  C  (|Dewll»  +  lDew2l-)2n  lWl"W2l2  M2  ’ 

where  we  have  bounded  the  Hessian  (of  3  9  -*■  F(  |  0 1 2 )  9  t  r^)  at  some 
intermediate  point:  C(|Dee|J2n  *  C(|D£W1|is+|De  w2|aj)2n  .  Combining  (11)  and. 
(12) ,  we  get 

(aVAvn)(VV 

( Au-Av  ) (u  -v  ) 

N  N  N 

c(  lDeul.  *  IV»I.)2"Iu-vnI2Ivv»I2  • 

which  with  the  triangle  inequality  establishes  the  following  estimate 

^3)  I u— un 1 2  *  C(!DeuL*  lDeVN^®)  I u_ VN 1 2  * 

for  any  vN  £  VN  H  W1  ’ *  provided  u  £  W1  ’ 2q  j”2rt  W1  ’ °°  * 

We  hereby  close  our  remarks  about  properties  of  Galerkin  solutions  and 

turn  our  attention  to  the  selection  of  the  basis-  or  Ansatz  functions. 


1  1 


4.  Asymptotic  expansion  (for  the  selection  of  basis  functions). 

Let  uN  be  the  Galerkin  solution  in  (given  by  (8))  of  the  variational 

problem  in  (9).  In  the  next  section,  we  shall  select  the  basis  -  or  Ansatz 

functions  { ip . }  such  as  to  establish  an  optimal  order  of  convergence  of 
J  j=0 

uN  to  u  (both  dependent  upon  e  )  as  e  tends  to  zero.  To  this  end  we 
employ  in  the  present  section  the  method  of  asymptotic  expansion  and  formal 
matching  to  "solve"  the  scaled  version  of  (4).  That  is,  expand  u  as  follows 

oo 

(HO  uU.n)  =  I  u(i,j)u,n)  ep(i,j) 

i  .j  =0 

for  {u(i'j)}.  .  ,  and  p  to  be  determined.  These  results  will  then  be 

i  »J 

applied  in  the  next  section  in  settling  said  convergence  as  e  tends  to 
zero.  It  will  turn  out  that  the  basis  functions  to  be  selected  are  the  ones 
predicted  (from  the  form  of  u^1*^)  by  the  asymptotic  expansions.  Recall  that 
formally,  (9)  and  (4)  are  equivalent  by  an  application  of  the  divergence 
theorem. 

Lemma  1 .  The  initial  exponent  and  term  in  the  expansion  (14)  are: 

(i)  For  u  <  0  : 

p(0,0)  =  -y  ,  u^’^i^.n)  =  CQU)  satisfying 

Co  '  "B  ’  Co(0)  =  Co(1 }  *  0  1 

(  i  i  )  For  y  *  0 : 

p(0,0)  =  0  ,  u^0,0\?,n)  «  CQ(£)  satisfying 

c0"  +  ((CQVn*V  =  -6  ,  CQ(0)  =  Co(D  -  0  ; 

1  2 


p(°,0)  =  -  2^-  .  u^0,0^U,n)  =  CQ(C)  satisfying 


((C  ')2n+1 ) '  =  -8  ,  C  (0)  =  C  (1 )  =  0 

O  0  0 


Also,  the  first  rpdependent  term  has  exponent  2  +  p(0,0)  in  all  three  cases 


(i)  -  ( iii ) 


Proof :  Let  v  =  u(0»0)ea  (a  =  p(0,0)).  Substitute  v  for  u  in  the 
scaled  version  of  (4): 

°£  (1  ♦  | D£v | 2n)  D£v 

.  u(0,0)  £a-2  +  (u(0>0)  ^n+1  (2n*1)(a-1)-1 

nn  n  n 

0  +  higher  order  terms 

with  boundary  conditions  at  n  *  ±  1 : 

u(°,0)  »-1  .  (u(0.0),2n.1  £(2n+1 ) (a— 1 )  .  h.0.,„ 

n  n 

Distinguishing  between  whether  a<1 ,  a=1  ,  or  a>1 ,  one  learns  in  either  case 
from  the  differential  equation  that  does  not  vary  with  n  ,  and  then 

from  the  bounda  -y  conditions,  u^^^  =  0,  such  that  u^,0^(5,n)  =  CQ(S)  • 

To  conclude  the  remainder,  let  the  first  rrdependent  term  be  called  u  and 

set  v  =  CQ  ea  +  u^’J^eb,  p(0,0)  =  a  <  b  *  p(i,j)  and  substitute  once  more  for 


D  (1  ♦  (C')2e2a  +  2C')u(iIJ)  ea+b  ♦  (u( l,J } )2e2b+(u( 1 ’ j  } )2e2b"2)n) 
to  o  t,  E,  n 


(C  V  .  u(1:J)  Eb  ,  e6”) 

o  e  n 


0  +  h .  o .  t . 


with  the  boundary  conditions  at  n  »  ±  1 

(1  +  ID  v|2n)  u(i,j)  =  ±  0(0  e2_b~w  +  h.o.t. 

1  e  1  n 

Distinguishing  between  whether  a  <  0,  a  »  0,  or  a  >  0,  one  learns  after 
lengthy  calculations  (formal  matching)  in  each  case  that: 

if  a  <  0:  then  y  =  -  (2n+1)a  >0,  b  -  2  +  a,  and 

((Co’)2n  +  V  -  -6  with  CQ(0)  »  CQ(1)  -  0; 

if  a  =  0:  then  y  »  0,  b  -  2,  and 

co"  +  ((  Co')2n+1)'  -  -0  with  CQ(0)  -  C0(1)  -  0  ; 
if  a  >  0:  then  y  -  -a  <  0,  b  «  2  ♦  a,  and 
C  "  -  -6  with  C  (0)  -  C  (1 )  -  0  . 

O  0  0 

Now,  distinguishing  between  a  <  0,  a  -  0,  and  a  >  0  in  each  of  the 
cases  (i),  (ii),  and  (iii),  the  proof  is  completed. 

//// 

Theorem  2.  Let  y  <  0  and  0  satisfy 

0  (  C°°(0,1)  and  V  k  f  N,  0(2k)(O)  =  0(2k)(1)  «  0. 

Then- the  exponents  in  (IH)  of  the  asymptotic  expansion  for  u  are  given  by 

05)  p(i  i  j )  ■  2i-(2j  ♦  1)  y  ,  yi.j  £  N; 

the  basis  functions  are  polynomials  in  n  of  even  degree  for  fixed  £  ,  and 

are  determined  by  the  recursive  formulae  (17)  and  (18). 


Proof ;  Define  [i,j]  -  21  -  2j \i  -  p(i,j)  +  u  ,  i,j  (  N;  we  concatenate 
terms  in  (14)  according  to  the  following  rule  (equivalence  relation): 


(i1  .  J1 )  (i2.J2) 

u  -  u 


<->  [i, ,j1 ]  -  [i2,j2] 


*1  =  i 2  Ji  *  J 2 


11  ~  i2 
J2 


u  if  J,  *  J2 


If  we  differentiate  through  in  the  scaled  version  of  (4),  we  get: 


[2F*( T)Up  +  F( T) ]u  +  [  ^  F'(T)u_u  ]  u_ 

5  SS  £2  t,  n  Cn 

+  C  K  F'(T)u2  ♦  rr  F(T)  3  U 
e4  n  e2  nn 


where  T  denotes  |Dru|2  *  u2F  +  e_2u2  and  F  is  defined  in  (2). 

^  n 

Define 


,(i  *  j) 


,u.j) 


[i1+i2.j1+J2]  *  Ci »J ] 


Ci1  ,J1  )  ( i 2 » J  2 ) 

u5  ^ 


l  u  (11 }  u  (i2,J2} 

Ci1+i2.J1+J2]  *  Ci.J]  5  n 

C i2  *  J  2 ]  *  Cl , 0] 


,(i » j ) 


[i1+i2.j1+J2]  *  Ci.J] 

[i1  ,j1  ]  S  2  ,  [i2,J2]  2  2 


(i1,j1)  ( i  2 » j  2 ) 


With  these  we  express 


(u  )2  =  e~2wN  .  where  N  =  I  v  (i,J) 
*  i»j*“0 


-2u  2 

u.u  =  t  e  0  ,  where 

S  n 


l  0(i*J)eCi*j]  -  e2  l  o( i+1 ' ^ >  eCi > j ] 


Ci » j ]£2 


[i+1 » j ]-2 


-  .2  -2p  4., 

(u  )  =  e  e  N-,  where 

n  2 


eV-  l  v2(’J>  eCi'J]  -  e4  I 

[l,j]*4  2  [i+2,j]£4 


„  (i+2,j)  [i.j] 
2 


Substituting  these  in  0  6),  we  obtain 


■(2n+1  )u  j  [2n(N  +  G2N  )n  1  N  +  ( N.  +  e2N  )n]  x  [  £  u^,j)  e[i,j]] 

'  2  1  i ,  j=0  ^ 


+  [4n(n  +e2N?)n  ^o]  x  [e^  l  u^  ',J>  e 
1  2  [ i+1 ,  j  ]£2 


2  l  u(i/1J)  e[i'j]] 


[2n(N  +e2N  )n*1N2+e"2(N  +e2N2)n]  x  [e2  £  uJ*+1,J)  e C 1  ,J  ]  ] } 

122  12  [i+1,j]*2  nn 


♦  c““  {  [  I  e[i’j]] 

i  .  j=0  w 


+  [  i  eCi,J]]  | 


[i+1 i j  J^2 


where  we  have  used  Lemma  1.  If  we  multiply  by  e*‘2n  and  list  terms  in 


n'  ”ji  kxilm  *  j.  - k  ^  w  *> 


descending  order  in  terms  of  upper  index  pair  for  a  given  power  [i,j],  we  get 


(17)  I  *  I  u(i.j*n)  £[i,j]  ,  „  ,  o  . 

[i+l.j-n]22  nn  Ci.j+n]*)  66 

R  is  an  abreviation  for  the  first  {...}  in  the  previous  equation  and 

contains  as  leading  term  (C  ')^n  u^i+1’^  .  if  we  multiply  the 

o  nn 

boundary  conditions  through  by  £(2n+1)y  1  ,  we  get 


(18) 


.  J 

.i+lfj 


[i+1rj+n]£2  n 


(i  +  l,j+n)  [i ,  j  ] 


+  [N.  +  e2N  ]n  l  u  (i+1,J)  e[i,j] 

1  2  [i+1.J]M  " 


=  ±3(5)  e2nW  at  n  =  ±  1 


Fix  [i,j].  Solving  the  Neumann  problem  (17)  and  (18)  for  u(i+1>j+n)  gives 
rise  to  a  solvability  condition  at  n  *  ±  1  yielding  an  ordinary  differen¬ 
tial  equation  in  5  fixing  an  additive  function  in  u^i,j+n^.  It  is  in  this 
O.D.E.  with  its  0  Dirichlet  boundary  conditions  that  the  conditions  on  B 
become  necessary  and  sufficient.  The  necessity  is  most  easily  seen  via  an 
even  extension  about  5=0  and  5=1.  Together  with  Lemma  1  this  deter¬ 
mines  u^'J)  uniquely,  i,j  N.  The  u^’^  are  even  polynomials  in  n  for 
fixed  5  ,  in  fact  induction  using  (17)  —  (18)  yields  that  u^’^(5,*)  is  of 
degree  2N  for  2NS[i , j ]<2N+2. 


5 

i 


i 

s 


//// 


We  remark  that  Thm.  2  would  hold  for  u  =  0  as  well  if  we  merely  drop 


the  j-term  in  p(i,j)  and  combine  the  leading  term  in  R  and  the  leading  term 

before  C  in  (17)  to  (1  +  C  2n)  u^i+1^  and  similarly  in  (18).  Now  for 

o  nn 

u  >  o. 

Theorem  3.  Let  \i  >  0  and  let  6  satisfy 

8  £  C°°(0,1)  and  V  k  t  H,  6(2k)(0)  =  8(2k)(D  =  0  . 

Then  the  exponents  in  (14)  of  the  asymptotic  expansion  for  u  are  given  by 

(19)  p(i,j)  =  2i  +  u  ,  V  i,j  t  N  ; 

the  basis  functions  are  even  polynomials  in  n  for  fixed  5  »  and  are 

determined  by  the  recursive  formulae  (20)  and  (21). 

Proof;  The  same  as  the  proof  of  Thm.  2,  except  that  the  leading  term  in  (17) 

now  is  the  one  from  R  ,  *.C  )2n  u^i+1’^  and  the  two  first  terms  in  (17) 

o  nn 

now  have  indices  (i+l.j-n)  and  (i,j-n)  respectively. 


(20) 


[i ♦1 t j ]-2 


l  (C  'l2"  .  I  (C  ')2n  cC1,j:l 

*•  <i  nn  O  55 


o  nn 


i » j  =0 


♦  R  .  I  *  I 

r«.ia _ n-.obn  r  <  <  . 


[i+1 , j -n]22 


[i,j-n]20 


(i,j-n)  [i  i  j  ] 

u55 


0  , 


where  R  is  R  excluding  the  two  first  terms  of  (14)  with  p(i,j)  given  by 


(19).  Similarly,  we  get  for  the  boundary  conditions  at  n  *  ±  1  : 


(21  ) 


[N^e^N^]0  I  u  :w  e 


(i+1 , j)  „[i,j] 
[i+1,j]>2  n 


l  u(i-1 ,J~n)  £Ci,j] 

[i+1,j-n]>2  n 


±  6(5)  • 


Again  we  have  to  concatenate  the  series  according  to  an  equivalence  relation: 


(VV 


(i2,J2) 

u  <=> 


<=> 


C i 1  » j 1 U  =  Ci2,J2] 


if  J, 


Also  here,  u^ (5,  •)  is  a  polynomial  in  n  of  degree  2N  for 
2N  S  [i ,  j]  <  2N  ♦  2  . 

//// 

We  remark  that  the  above  method  is  quite  general.  It  could  for  instance 
be  generalized  to  a  layered  material  where  F  might  depend  also  on  n  .  In 
such  cases  we  would  no  longer  get  even  polynomials  in  n  as  Ansatz  functions, 
but  the  same  type  of  recursion  formulae  would  determine  u^’-^'s  n-dependence 
and  the  p(i,j)  would  remain  the  same. 


5.  The  selection  of  basis  functions. 


We  will  now  apply  the  results  from  the  last  section  to  the  selection  of 
the  basis-  or  Ansatz  functions  }J_q  •  This  selection  fixes  the  sub¬ 

spaces  VN  as  defined  in  (8),  and  we  let  u  and  uN  be  the  solutions  of  (7) 
and  (9)  respectively.  The  aim  of  this  section  is  to  make  said  selection  so  as 
to  have  optimal  order  of  convergence  of  u-u^  as  e  tends  to  zero  at  fixed 
N  . 

Before  we  state  and  prove  convergence  of  u^  to  u  for  e  tending  to 

zero  in  the  three  asymptotic  ranges  characterized  by  the  sign  of  p  ,  we  let 

2  2 

D1  be  the  operator  defined  by  D^u  =  (d  u/d£  )  mapping  Dom(D.|)  = 

w2,2n+2(o,i )  n  W1,2J+2  (0,1)  ->  L2n+2(0,1)  . 

Theorem  4  Let  p  =  0  and  n  £  Z+  . 

\/  N  6  H,  3  CN  independent  of  e  such  that,  if  8  6  DCD^)  then 


2N+2  - 


2n+2 


IUUnI  ^  CN  e 


(2n+1 )‘ 


(n+1 ) (2n+1 ) 


a 


Proof:  For  convenience,  let  n  =  1.  Let 


N 


where 


Jo  “J"'2' 


UjNU,n)  =  ^i=0  XiN  (5)*i(n)-  171611 


D  UM 

e  N 


r  v  2 1  r  21-1 

l  l  (uiN)E  e  .  I  (u  )  e  J 

j  =0  jN  5  j  =0  jN  n 


20 


m 


I,*J 


iMI 


■•ft! 


Is  'V 


\»,y 

M 


i 


SVC 


<",*1 


i  I  f  J, 

V’’.k 


i 


V.y 

I 


1 


m 

11 

« 


»»;a 

»«*>l 

•<? 

I?3 


s. 


m 


t.\v 


i  4  14  14 

On  the  grounds  of  denseness  of  W  ^  (0,1)  8  W  ’  (-1,1)  in  W  ^  ( Q)  and 


—  /  1  4  * 

the  fact  that  Au„  t  (W,!.(Q))  ,  we  restrict  v  to  be  of  the  form  c  •  d. 

N  (0) 


where  c  t  W^^OJ)  and  d  (  W^’^-I.l)  ,  for  convenience.  Then  (9) 


becomes  with  u. 


I  la  ^  lDe  unI  ^  J0  (ujN^  e  J  ’  J0  (ujN)n  £  J  ^ 


x  (cd,  e  cd')  d£dn 


=  /  6(0  cd  d£ 

r. 


=  2  J  6  cd£  d(1 )  , 
0 


for  any  such  v  which  is  even  in  n  (we  shall  let  v 


u-  u^  ,  even  in  n) 


and  where 


’  *  iv/ 


““o/n3  E'2 


^  +  2(u0N)  (uIN)  -1 

5  n  n 


l-2(uON)  (uiNJ  +  2(uON)  (U2N}  +(U1N}  ^ 

5  5  n  n  n 


2  i  4N-2 


^2(UN-1  .N^.^NN^  +  (UNNJ  ^  E 

5  5  n 


r ,  .  2  4  N 1 

+  UNN  ^  E  ]  • 


V 


asymptotic  expansion  (possible  if  8  6  C<11'<(0,1)  with  6i'£:k;(0)  =  6  (1)  =  0 

for  OSkSN)  ,  see  Thm.  2  with  subsequent  remark.  However,  the  coefficient 

..  .  2N  +  1 

equation  for  e  : 


S  !  2)  UuNN)  c*d)  +  .  .  . 


(the  coeff.  to  e2N  in  Id  u.,|2)[(u  )  c'd  +  (u  )  cd')d£dn 

1  e  N‘  ON  _  IN 

C  n 


up  to  the  coefficient  equation  for  e°  (highest  power  present): 


/  /  3  c’d  di*dn  =  0 


fl  NN  C 


will  change  as  we  include  more  basis  functions.  Define  u^  (modulo  an 

N  2 

additive  function  of  £  )  to  be  that  exact  solution  (  l  u  t  J)  of  the 

j=0 


-3  1 

coefficient  equations  of  e  ,  e  , 


f  1  N 

{u.,,}  is  thus 
JN  j-o 


/  N_  I 

defined  iteratively  provided  8  K  D(D1  )  ,  see  Thm.  2.  Furthermore,  if  we 

suppose  8  6  D(D.|N)  ,  the  coefficient  of  e2N  +  1  in  Au^tv)  can  be  written 


/  la  "  ♦  <uo»)?2)(<uN»i,ii*i)n  od'> dMn 


i  y  i  -  y 

where  we  write  cd'  as  (  —  v  )e“  e 

t  n 


Thus,  V  v  6  (£J) 


(Au  -  Au„) ( v)  =  G(v)  -  Au„( v ) 


A 

$ 

$ 

£ 


G(v)  -  (G(v)  +  R(un  ,  v))  , 


where  (for  e  bounded;  I R I  -*•  ®  as  c  for  e  -*■  °°) 


n/- —  v  |  -  ii  2N+2-1/ 

R<V  v)l  «  CN  lvl  e 


which  together  with  strong  monotonicity  implies 


a  lu'uN|4  *  |(Au-AuN)(u-uN) 


.  _  ,  -  ,  2N+2-V 

s  cnI  »l  E 


With  (10)  this  implies  (with  a  new  constant  CM) 


lU-“Kl  S  CN  C 


2N+2-l/u  1 

9  +  6 


for  n=1 .  Larger  n  change  only  the  trivial  and  tiresome  algebra. 


We  immediately  obtain  the  following  estimate  in  W 


Corollary  1  Under  the  same  hypotheses  as  in  Theorem  3  and  supposing  that  u, 


u  W 
N 


independent  of  c  such  that 


mm 


Iu  unI2  s  cn  c 


2N+3/2 


ft 

& 

!: 


•i* 

»« 

I 

| 


% 

i 


I 

I 

| 


I 


1 


Proof :  Again,  let  n=1 .  Now  string 


Iu-UNI22  *  (Au-Aun)(u-un) 


*  |R  (uN,  u-uN)|  i 


R(un,v) 


U  In  ^  +  (uON} (uN+1 , N+1 }  cd’  d5dn 


+  e4Ni  I  /Q  (UNN) ^  C’d  d5dn 


5  CN  ^  lvl2  ’  V  v=Cd  *  w(0)  • 


inequality  (13)  for  v„  =  u  together  to  obtain  the  estimate;  w.l.o.g.  we 


have  let  v=c  *d  ,  c  W1,^q+2(0,1)  ,  d  f  W^’2n  2 ( — 1  *  1 )  . 


The  next  Thm.  covers  the  case  u<0  . 


Theorem  5.  Let  p  (  R  n  f  Z  +  .  V  N  (  R,  3  C„  independent  of  e  such 

that  if  8  (  DtD^)  then 


e(N,  pi) 


+  d(n) 


!u-un1  s  cn  e 


(2n+1 )‘ 


where 


e(N,y)  -  2N+2-U- 


2n+2  ’ 


ft&V/V'-Y' 


m 


i.«*’ 


i 


Si 
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IS 
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m 
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mmmm 
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(n+1 ) (2n+1 ) 


Proof :  Let  (  [i,j]  »  2i  -  2ju  ,  i»j  *  M) 


[i , j ]<2N+2 
i  *  j  €  N 


..  (i.j)  r  Ci,j]-W 

N  E 


where  u^’^U.n)  =  x^1  (5)'l'k(n)  ,  where  (in  case  u  £  Q)  we  take 

the  representative  (i.j)  for  which  -ju<1  .  The  remainder  of  the  proof  is 
then  a  mere  imitation  of  the  procedure  shown  in  the  previous  proof  with  the 
aid  (as  also  there)  of  Thm.  2. 


For  u>0  ,  we  need  an  extra  assumption  on  B  : 


Let  B  (  (W1,2n+2)*  »  W  2n+1  satisfy: 


The  function  CQ  £  W1,2q+2  (0,1)  defined  a3  the  solution  of 


J  {(C^)2n+1c'  +  Bc}d£  -  0,  Vc(m1,2J+2  (0,1)  , 


(A)  satisfies  that  there  exists  u  ^’^o)  ^  *  solution  of 


I  I  (c’0)2nVn  =  g(v)  ’  Vv  ^w1,2(0)  (fl) 


.  1  2n+2  * 

for  any  given  g  k  (W  ’  ^  (55)) 


Theorem  6  Let  p  €  R  ,  n  £  Z  .  Let  8  satisfy  assumption  (A).  V  N  £  N, 


/  N  -3 

let  8  t  D(D. )  .  Then  3  CN  independent  of  e  such  that 


I I  ^  e 


e(N,y)+d(n,y) 

/  . 4  \  2 


(2n+1 ) 


where,  again 


e(N,u)  =  2N+2-y  -  2^2  '  and 


d  (n,  y)  =  2n( 1 -y  - 


Proof :  Let  (with  the  notation  of  Thm.  3) 


,a.j) 


[i,j]  - 


[ i , j ]<2N+2 
i,J  *  N 


(i.J) 


,(i  *  j ) 


tj)k  ,  where  (for  y  Q)  we  take  that 


where  u^1’^  *  Ik*Q  xki,d  \  ’  where  (for  y  Q)  we  take  that 

representative  (i,j)  for  which  2n+T  w  <  1  ♦  Otherwise  the  proof  is  like 


that  of  Thm.  4  with  the  aid  of  Thm.  3  and  the  inequality  (10)  for  y>0  . 


Assumption  (A)  is  rather  restrictive  and  can  be  dispensed  with  through 


containing  ourselves  to  local  estimates  or  estimates  in  weighted  norms. 


As  presented  in  Corollary  1  for  y=0,  it  is  possible  also  for  y*0  to 


12  1  2n+2 

obtain  estimates  in  the  -  norm  rather  than  the  w  ’  (q)  “  norm> 


Corollary  2  Let  y<0.  Under  the  same  hypotheses  as  in  Thm.  4  and  supposing 


✓*  1  ® 

that  u  and  u^  (.  W  ’  ,  3  C^  independent  of  e  such  that 


"V  ?  : 

w 


$ 

wl» 

m 

SvV 


i 

l 


-  O.  -  r 


I 


where  e(N,u)  is  the  one  defined  in  Theorem  4  with  n=»0. 
The  proof  is  like  that  of  Corollary  1. 


Corollary  3.  Let  p>0.  Under  the  same  hypotheses  as  in  Thm.  5  and  supposing 
_  u  _ 

1  ,« 


«.  .  2n+1  2n+1  r 

that  c  u  and  e  u..  c  W 

N 


C^  independent  of  e  such  that 


et»,u)  -  f^r 

s  cN  « 


where  e(N,u)  is  the  one  defined  in  Theorem  5  for  n=0  . 
Proof ;  As  in  Corollary  1  , 


Iu-uj2  S  C„  £ 


e(N,p) 


but  (13)  now  states 


|u-u„|  S  C  <|Du|.  •  |IWN|  )2n  |u-uN| 

2  *2 

2ny 

r  2n+1  I  -  | 

S  C  e  |u-uN| 


//// 


The  combination  of  the  two  inequalities  completes  the  proof. 

It  1 3  important  to  note  that  we  have  made  the  selection  of  basis  or 

r  i  N  — 

Ansatz  functions  14^}  ia0  (and  thus  V^)  by  constructing  u^  out  of 

^UN  ^ [ i » J  "•  *  w^ere  UM  ^  (5»h)  *  I  u_n  x„  (5)  (h)  are 


I  ] <2N+2 

determined  through  solving  equations  which  are  the  weak  variants  of  the  formal 


matching  equations  already  encountered  in  Section  4. 

We  observe  that  for  our  model  problem  (the  particular  form  of  F  given 

i  i  N 

in  (2)),  we  obtain  the  same  selection  of  (namely  the  even  poly¬ 

nomials)  in  all  three  asymptotic  ranges  ( p=0  meaning  finite  limit  load, 
u<0  :  zero  limit  load,  and  u>0  :  infinite  limit  load).  This  would  hold 

for  more  general  F  (see  the  paragraph  following  (9))  including  such  that 
have  a  factor  depending  upon  n  (e.g.,  layered  materials);  here  the  selection 
of  would  be  made  according  to  the  same  scheme,  only,  the  sequence  of 

second  order  differential  equations  determining  would  have  coefficients 

depending  upon  n  and  would  no  longer  merely  be  polynomials. 

In  the  next  section  we  intend  to  show  that  this  choice  is  optimal  with 
1  2 

respect  to  "  norm  convergence  under  appropriate  boundedness  conditions 

on  D  u  ,  D  u„  in  L°°  . 
e  e  N 
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6.  The  optimality  of  above  selection  of  basis  functions  (at  fixed  N) ■ 
Define  for  *  W1  ,2n+2(-1 , 1 )  , 

WN  “  {  E  Vi  :  ci  W1,2nQ2  (0,1)}  , 

I  i  N 

and  let  lVi«0  be  the  basi3  functions  selected  in  the  previous  section 

(defining  VN) .  There  we  established  the  following  convergence  rate  in 
1  2 

W  ^  -norm  for  e  tending  to  zero: 

IU_UnI2  *  CN  ¥(e) 

under  certain  assumptions  (see  Corollaries  1,2,  and  3),  where 


2N+|  -  u 


,  for  u£0 


2N+|-  -  u  _  2n+T  u  *  for  M>0 


r  i  N 

Regarding  the  question  of  optimality  of  the  selection  of  I41  f  we  now 

1  0 

state  and  prove 

Theorem  7.  Let  u  S  0.  Under  the  hypotheses  of  Corollaries  1  and  2,  suppose 


{$ .}N  exists  such  that  (for  all  N  (.  M) 

^  J  A 


inf  |u-w  |  S  $(e) 
WN  ^  WN  2 


for  some  function  4>  with  the  property  that 


is  bounded  for  all  e  >  0 


Then  WN  2  VM  . 


Proof :  The  proof  is  by  induction  on  N.  Take  the  particular  function 


vM  -  l  u  related  to  u..  of  Thm.  H  and  5. 

“  [ i , J ]S2N  N  N 

such  that 


lvN  "  «N|  *  *^e) 


'  wk-  e  WM 

N  N 


In  particular,  let  N  =  0: 


lv0  "  woL  =  K  6 


(0A0)  e-“  -  w, 


l  -  v 

w0|2  =  0(e?  ) 


We  thus  get 


|u(0’0)  -  wQ  =  0(e2) 


such  that  (if  we  denote  w^0^  =  wQ  e+M  ) 


(0,0)  _  (0,0)  2  .  ,  (0,0)  _  (0,0)  ,2  „  „  3 


dx  S  |u->-  -  I  *  c0 


whence  w(0^0)  =  u(0^0)  and  uQ  =  wQ  £  WQ  . 


Suppose  next  that  1  S  span  (<|)1}N  ,  then  there  exists 


WN  ^  WN  3uctl  that  is  defined  in  Thm.  2) 


I'n  ‘  "»l. 


*  I  1  “H 

[i,J]S2M 


u  (i.J)  Ci tj ]~u  _  -  j 

N  e  N1. 


0(e2N+f-  »*)  . 


where  um  *  2.  u  n  xu  K  as  be'ore*  This  means  that 

N  K  K  K 


I.  I  Un-J>  En.j]  -  » . 


[i,j]£2N 


\*‘rx  ^  t  j„  ■i,,0,vaH\1  i  “n1,j>eCI,j]'“-  “»>' 

k«0  L i » J  J<2N 


0  (e2N_U+2)  , 


r  1 N 

such  that  upon  taking  the  limit  e  *  0  ,  span  i ♦  ^ } i*o  *  This  completes 


the  proof. 


We  note  that  it  is  possible  to  trengthen  this  Thm.  to  allow  4>  for 


which 


eP  S  o(1)  ,  as  £  *  0  for  any  p  i  [0,2], 


by  employing 


e_1/2  CJJ  1 (vVl2^  *  Ivn_wnI, 


0(e2N-pTP) 


This  means  that  even  if  we  were  satisfied  with  a  convergence  rate  of 

2N+|  -  u  -  p  ,  p  £  [0,2]  ,  we  would  still  need  to  use  {i|>.  (N  .  We  shall 

2  i=0 

employ  this  technique  in  the  next  theorem. 


ifW 


Theorem  8.  Let  y  >  0.  Under  the  hypotheses  of  Corollary  3.  suppose 
exists  such  that  (for  all  N  (  M) 


inf  |u  -  wN|  £  4>( e ) 

WN  *  WN  2 


for  some  function  $  with  the  property  that 


Jjn+1 

y|yy  e  2n+1  =  o(1)  ,  as  e  +  0,  for  some  p  +  y  i  [0,2] 


Then  W„  V  . 

N  N 

Proof :  For  p  +  U  ~  [0,y  ]  ,  the  proof  would  merely  copy  that  of  the 

r  i  •  i  JJ 

—  (i  j )  ~  2n+1 

previous  Thm.  with  j]^2N  UN  *J  e  and 


[i.j]  =  2i  +  ^y  . 


As  is,  let  N  =  0 : 


i  §  (0,0)  2n+1  I  f 

IV  woL  =  >u  o  e  '  W0I  =  0  U 


-  u-  2^  y  -  P  -  2 


«n+l  ,  3 

r  "  2^T  y  '  P  +  2^ 

•  o  [t  J  , 


^  If  \h  K 


(0,0) 


2  ‘4 

-oil  ] 


Un+1  „  3 
_  y  + 


Since 


(0,0)  =  Q  ^  letting  e  +  0  ,  <j> _  *  0,  and  i|»n  £  span  Ur 


3n  0 


With  this  modification,  the  induction  step  will  carry  over  from  Thm.  7. 


I 


I 


i 


% 


i 


ft 

ft 

ft 


ft 

ft 

ft 


$ 
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This  ensures  optimality  for  n=0  and  otherwise  p  <  2(2n+1 )/(4n+1 )  . 

00 

We  emphasize  that  for  p  >  0  we  used  an  L  bound  on  the  gradients  of 

_  _P _  _  _P _ 

2n+1  2n+1 

e  u  and  e  u^  ,  which  implies  that  we  have  an  e  -dependent 

constant  CCjDul^,  JDvJ^)  in  (13). 

In  Thms.  4-6  we  established  convergence  rates  for  e  tending  to  zero 
with  N  fixed.  Even  though  we  are  dealing  with  "thin"  structures,  t 
nevertheless  has  some  finite  value  over  which  we  may  have  no  control.  Even 
if  e  is  "very  small"  we  have  seen  that  convergence  rates  are  limited  by  the 
degree  of  smoothnesss  of  the  load  6  .  It  is  therefore  of  interest  and  im¬ 

portance  to  prove  convergence  for  the  number  of  reduced  models  (N)  increas- 


7.  Optimal  convergence  for  Increasing  number  of  reduced  models. 

Let  e>0  be  fixed.  Let  uN  be  the  Galerkin  solution  in  VN  (given  in 
(8))  of  (9)  as  before.  Proving  convergence  of  u-uN  for  N  tending  towards 
infinity  reduces  to  an  approximation  theoretical  problem  because  of  (10).  Fo 

y 

completeness  we  outline  it  here,  cf.  [31],  Babuska,  Szabo,  and  Katz  [3],  and 
Canuto  and  Quarteroni  [5]  for  L^-theory  and  Meinardus  [23]  and  Jackson  [ 1 3D 
for  L°°-theory.  We  then  interpolate  in  order  to  obtain  L^n+^-estimates. 
Details  are  found  in  [14], 

The  solution  u  is  an  even  function  in  n  so  we  define  for  any  p 


[2,  •]  ,  m  f  Z+ : 


-  m,p 


W  (-1,1)  =  Wm,P(-1,1)  {v:v  even} 

.(1  ,m)p  *1  ,p  „m,p  -tn.p 

Let  w(^  n)  *  W  (v  t  v^  t  }  ,  2£p£°°  ,  m  t  Z+  ,  where 

»  m.P 

v  C  W  if  v  and  its  first  m  partials  w.r.t.n  belong  to  Lp  in  the 


usual  generalized  sense.  Define 


L(1  ,m)p 

V.n) 


|v|  1  +  |am+1v/H9nm|  +|3mv/9nm| 


Let  PN  be  the  space  of  functions  which  are  polynomials  in  n  of  degree  at 
most  N  and  PN  be  those  functions  in  PN  which  are  even  in  n  • 

Using  fairly  standard  argument  along  the  lines  of  the  references  cited 
above,  one  obtains: 

Theorem  9.  For  all  N  £  N  ,  m  t  Z+  ,  there  exists  an  operator 


,0  ,m)2 
A  .  u 

2N  •  U.n)0 


PnM  -  w  (n)  £  w  (n) 

c  N 


V 

*5 


which  is  a  continuous  linear  mapping: 


r  .(1,tn)2  .1,2 

I  -A  eA  (W,c  nv  ,  W  (0))  ,  and 

2N  5 ,  n) 


.(l,m)®  .1,® 

1  -A  2N  ^  •  "  <!1))  • 


with  norms  bounded  by  a  constant  times  (2N)'  m  . 

We  then  obtain  the  main  result  of  this  section. 


Theorem  10.  Let  u  w(1*m)2n+2  (0)  for  an  m  f  Z.  ,  n  t  Z+  .  Then  there 
-  U,n) 

exists  a  positive  C  independent  of  u  and  N  such  that  for  N  t  Z+ 


l2fv  ^u~vly  l,2n*2  S  C(2N)'"’|u|^;J)  . 

V  t  N 

Proof .  Employ  interpolation  by  the  K-method,  see  Bergh  &  Lofstrom  [4]  Thm. 
6.45  (5),  we  obtain  from  Thm.  9  (choose  6  =  1/(n+1)) 


^  .(1,m)2n+2  .1,2n+2 

(i  -A  „)  **  <wUin)  ,  W  (n» 


with  a  bound  on  the  norm  proportional  to  (2N) 


This  yields  the  immediate 

Corollary  4.  Let  u^R»  n  f  Z+  .  V  N  f  Z+  ,  3  C  >  0  independent  of 

.(1 ,m)2n+2 

such  that,  if  uf  W  ,  c  .  (0), 

v  s  t  n  / 


|u-u,J  S  C  (2N)‘ 


u  being  the  exact  Galerkin  solution  of  (7). 


Proof:  Combine  (10)  with  Theorem  10. 


The  similar  Corollary  about  W1,2-rates  out  of  Thm.  9  is 
Corollary  5.  Let  p  (  R  ,  n  (  Z  ,  V  N  f  Z+  ,  3  C  >  0  independent  of 


~(1,m)2 

N  such  that,  if  u  (  K,.  ,  (n)  W 


1  00 

bounded  in  W  ’  uniformly, 


|u-uN|2  <  C(2N)  "m  . 

Proof :  Combine  (13)  with  Theorem  9. 

////  • 

We  will  be  in  the  position  to  show  that  the  subspace  VN  gives  quasi- 
optimal  rates  of  convergence  in  W1,2  (in  the  sense  of  n-widths)  after  the 
following  definitions  (i.e.,  in  general  one  could  only  improve  the  constant 
not  the  rate): 

The  (Kolmogorov)  N-width  of  W  (a  subset  of  W1’2)  in  W1  ’ 2 ( fj )  is  given 
by 

^..(WjW1’2)  =  inf  { E( W;  W..):W„  is  an  N-dimensional  subspace  of  W  } , 

N  N  N 

where 

E(W;  W  )  =  sup  inf  |x-y| 


x  £  W  y  £  W. 


is  the  deviation  of  W  from  . 


A  subspace  WN  of  W  of  dimension  at  most  N  for  which 


dN(W;  W  ’  )  =  E(W;  WN) 


i  2 

is  called  an  optimal  supspace  for  dN(W,  W  ’  ).  The  notation  follows  that  of 
the  monograph  by  Pinkus  [27]  . 

Let  W  =  Wm,2n+2( -1 , 1 )  .  In  said  monograph  [27],  the  N-width  is  com¬ 

puted  in  Chapter  VII.  Set  B  =  (v  £  W:|v|£lf  ,  then  there  exist  positive 
constants  C  and  D  independent  of  N  for  which 

D(2N)"m  £  d.( B ;  L2)  £  C(2N)"m  . 

N 

This  yields  immediately 

D(2N)1_m  £  dM(B,  W1 ,2(-1 ,1 ))  £  C(2N)1_m  . 

N 

In  view  of  the  denseness  of  W1,2q+2(0,1)  ®  W1  ’ 2n+2( -1  , 1  )  in  W1,2rlQ2(fi)  , 

this  shows  that  the  estimate  in  Corollary  5  is  optimal  (modulo  the  constant  C) 

and  that  this  is  obtained  by  the  subspace  VN  .  Thus  [i|>.  \  (or  any 

1  i  =0 

{<j>.  }N  with  the  same  span)  provides  a  quasioptimal  choice  with  regard  to 
1  i  =0 

convergence  at  fixed  e  for  the  number  of  reduced  models  (N)  increasing. 


8.  Computational  results 


Consider  the  numerical  problem  of  finding  the  dimensionally  reduced  solu¬ 
tion  of  (7)  within  V^.  Let  n=1  in  (2)  as  an  example.  We  will  consider  two 

2  2  2 

loads  8  :  8=1  and  8(x)  =  Air  (1+3A  cos  Trx)simix.  We  limit  ourselves  in 

advance  to  N  =0,  1,  and  2. 

2  4 

Select  4>q  =  1  ,  <J>2=n  ,  and  4>^=n  .  These  form  the  same  span  as  the  first 

three  <Jj  .  as  could  be  found  from  Thms.  1  and  2.  Let 


I  v0, (£)♦,,(  n) 


j=0 


2  j  2  j 


with  v2j  £  W1 ,8,^(0, 1 )  .  Substituting  this  expression  into  (7)  yields  a  system 

of  N  nonlinear,  second  order  ordinary  differential  equations  in  £  ,  upon 

f  i  N 

integration  w.r.t.n  over  [-1,1],  to  be  solved  for  iv2j‘j=o  •  We  convert 
this  into  a  first  order  system.  Let  w^+1  =  v2j  ,  and  wj+jj+2  =  v2j  for 
0£j£N  .  We  arrive  at 

A(w)w'  =  R(w, 8)  , 


.  ,  .  ...2N+2  ,  /  „2N+2x2N+2  .  /  „2N+2  ,, 

where  =  (w(j))j=1  ,  A  t  R  ,  and  R  c  IT  .  Upon  inversion 

of  A  ,  the  system  is  in  a  state  such  that  it  is  numerically  solvable  using 
the  nonlinear  ODE-system's  solver  written  by  V.  Majer,  cf.  Babuska  and  Majer 
[21],  and  Majer  [22]. 

From  that  point  on  we  have  taken  mainly  two  approaches.  The  first  to  be 
mentioned  is  applicable  in  general  for  nonsmooth  or  noncompatible  boundary 
loads  8  as  well  as  for  all  three  asymptotic  ranges  (y<0,  y=0,  and  y>0  ). 
The  second  is  applicable  only  for  smooth,  compatible  8  and  for  u£0  ,  in 
which  cases  an  asymptotic  expansion  can  be  generated  automatically. 


We  now  describe  the  performance  of  the  first  approach.  Let  S»1  (non¬ 


compatible)  as  an  example. 

The  exact  energy  (and  solution)  being  unavailable,  it  is  rather  found  by 

extrapolation  from  the  energies  E(uq),  E(u1),  and  E(u2)- 

It  is  intuitively  obvious  that  in  the  range  | D^u j  << 1 ,  E(uN)-E(u)  can  be 

interpreted  as  the  square  of  the  difference  in  the  f^-norm.  Also,  there  is 

1  to 

some  relation  to  the  difference  in  the  W  ’  -norm.  Define 

irl 


$:R2  )  p  -  $(p)  =  d4|p|2)|p|2  .  Then 


1 

E(un)-E(u)  =  2  J7  /  dp4>(Deu+s(D£.(uN-u)))  *De(uN-u)dsdCdn  -  G(un~u) 


fl  0 


and  since 


/  /  e  d  *(D  u)*D  (u-u)  -  G(u  -u)  -  0 

'  ‘  P£  eN  N 


(u  is  the  minimizer)  we  get 


E( u  ) -E( u)  -  / /  eD  (u  -u)H  D  (u  -u)  , 

N  2  '  e  N  «  e  N 


1  s 


where  H  =  [  \  d2  <J»(D  u  +  tD  (u-u))  dt  ds  with 

0  0  P  e  e  N 

2  2 

(d  4>).,  =  6 .  .  ( 1  +  I p I  )  +2p,p.  .  Thqs  we  can  estimate  with 

P  ij  ij  11  1  J 


1.2 


7*/toT 

2%— 


>  0 


and 


E(uN)  -E(u)  2  4  V  /  /  J  |De(uN-u)|2  ♦  1  |D£(un-u)|4 

E(un)  -  E(u)  S  c  J  /  j  |De(uN-u)|2 

+  3X2  {|Deuj2|D£(uN-u)|2  «■  I Dg(u^-u) j  4 | 


too 


m 


VJVJ 


If  we  furthermore  assume  that  uN,u  t  W1'®  with  ,  |D£;u|oo  <  M  for 
all  N  for  some  positive  constant  M  ,  then  said  relationships  between  the 
difference  in  energies  and  the  norms  of  the  energies  emerge: 


E(un)  -  ECU)  -  |Vu| 

H 


MVU!4  *  E(uN)  -  E(u)  S  (  \  ♦  15X2M)2|a|l/2|uN-u| 


The  latter  assumption  is  dubious  however  in  the  limit  e  +  0  for  u  >  0  (see 
Lemma  1).  With  said  assumption,  Kosler  [19]  proved  the  usual  Cea’s  Lemma  type 

estimate  | u-u  |  S  C  inf  ,  |u-v  |  i  •  See  03).  For  u  *  0  , 

H  N  t  N  N  H 

we  thus  expect  "linear"  behavior  (as  well  as  for  u  «  0,  if  M  <<1).  See 
Corollaries  1  and  2. 

The  expected  rates  of  convergence  (for  e  +  0)  measured  w.r.t. 
difference  in  energy  are  (cf.  Vogelius  &  Babuska  [30]  for  u-1 ,  linear  case) 


Linear 


(N-0,1 ,2) 


Nonlinear 


(N-0,1 ,2) 


u  =  -1 


(5,6,6) 


(3,4,H) 


0,2,2) 
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vg »  ^  ’  3 ' 


where  the  *  indicates  that  they  are  based  upon  W1 ’ ^-estimates .  If  we  used 

the  w1 ,2-estimates,  we  would  get  the  rates  for  the  linear  case  except  that 

12  2 

for  y-1  ,  we  get  (  -  z,  rr,  ■=■)  instead. 


u 

e 

E(uq)-E(u) 

E(u1 )-E(u) 

E(u2)-E(u) 

.1 

6.35  e-8 

2.93  e-10 

1  .35  e-12 

-1 

.01 

6.64  e-1 3 

2.80  e-16 

1.18  e-19 

.001 

6.66  e-18 

0 

0 

.0001 

7.00  e-23 

0 

0 

.1 

6.31  e-6 

2.92  e-8 

1.34  e-10 

0 

O 

1 

CO 

r\> 

2.07  e-7 

2.92  e-10 

4.12  e-13 

.01 

6.62  e-9 

3.90  e-12 

2.29  e-15 

10-5/2 

2.10  e-10 

0 

0 

1  . 

5.90  e-9 

6.44  e-10 

7.03  e-1 1 

+1 

.1 

9.92  e-10 

4.58  e-12 

2.11  e-1 4 

.01 

1 .04  e-10 

4.40  e-1 4 

1.87  e-17 

.001 

1.03  e-1 1 

1.00  e-1 4 

9.69  e-18 

Table  1 


Differences  in  energy  dependent  upon  e  and  u  for  0  =1  • 


The  last  column  E(u2)-E(u)  has  an  extra  factor  e  built  in  on  account  of  the 
extrapolation.  We  notice  the  expected  linear  rates  for  u<0  ,  to  some  extent 
u«0  ,  and  in  an  e-range  for  u>0  ,  whereas  there  is  an  expected  marked 
slowdown  for  y>0  as  e  +  0  . 

The  convergence  rates  for  N+®  exhibit  the  same  type  of  behavior.  From 
v  -4 

Vogelius  and  Babuska  [31 L  we  expect  E(uN)-E(u)  -  N  ,  see  Figs.  1,  2,  and  3. 
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The  higher  order  (N)  models  become  necessary  not  only  for  e  not 
small,  but  also  when  sufficient  smoothness  is  absent  (boundary-  and  interior 
layers).  Pertaining  to  the  latter  case,  we  also  need  to  show  extra  considera¬ 
tion  when  we  want  to  resolve  the  (from  a  practical  viewpoint  very  interesting) 
"nonsmoothness"  in  the  solution  well.  This  would  namely  mean  that  we  would 
like  to  resolve  w j ,  j  large,  well.  But  Wj’s  support  lies  in  a  (still) 
narrow(er)  band  near  the  singularity  (for  j/e  large  enough)  calling  for 
separate  adaptively  generated  meshes  for  those  Wj  .  For  the  linear  case,  see 
Vogellus  and  Babuska  [32].  See  Figs.  5,  6,  7,  8,  and  9. 


We  now  describe  the  performance  of  a  second  approach  applicable  for 

smooth,  compatible  6  and  y£0  .  In  this  case,  it  is  possible  to  generate  an 

asymptotic  series  for  the  exact  solution  along  the  lines  of  Section  4.  Taking 

a  sufficiently  large  number  of  terms  in  this  series,  it  is  thus  feasible  to 

make  direct  error  computations  in  the  relevant  norms.  As  an  example,  take 
2  2  2  2 

6  =  Att  (1 +3A  it  cos  irx)sinirx,  ,  where  A  is  a  continuation  parameter. 

(This  choice  corresponds  to  taking  CQ  =  Asirnrx,  cf.  Lemma  1  part  ( i i ) . ) 

Let  the  exact  solution  be  given  by  (cf.  Thm.  2  with  subsequent  remark) 

u  -  l  u(j)  e2j 
j=0 

where  each 

u(J)  »  1  a  U)n2k  *  J  £  N  . 

k-0  J’ 

Denote  by  uM  the  Mth  partial  sum  in  the  above  asymptotic  series.  Then 


I  KunI  -  Ivun|  I  s  |uu*| 


where  |U_UNI  is  sought  for,  Ium-unI  is  comPutable>  and,  by  inspection 


HUi.; 


o(c2M  *  3/2> 


and 


ou2M  *  5/,>  . 


From  Cor.  1  and  Thm.  3  it  follows  that 


and 


lU'Unll,2 

.4 


OU2"  *  3/2) 


2N  +  7/4  1 

ou  9  * 


To  investigate  the  asymptotic  rate  of  convergence  for  e  +  0  for  N-=0  and  1 
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respectively  it  will  thus  suffice  to  take  M»1  and  2  (for  the  W  ’  -norm)  and 
M=0  (for  the  W1,l4-norm). 

The  terms  u^  ,  j-0,  . ..,  M  ,  are  generated  via  SMP  (the  symbolic 
manipulation  program  by  Inference  Corporation)  using  the  recursive  formulae  in 

the  proof  of  Thm.  2  with  the  twist  that  p=0  .  |uf-u„|  is  then  computed  via 

1  M  N 1 

Simpson's  rule  over  51  points  in  [0,  1/2]. 

On  the  next  two  figures  (Figs.  10  and  11)  are  shown  the  errors  in  the  two 

names  W1'^  and  W1’1*  for  N=0.  We  see  the  optimal  rate  for  lu-u.,1.  _  and 

1  N 1 1  f  d. 

observe  a  rate  for  lu-u^^  ^  ,  which,  not  surprisingly,  is  superior  to  the 
rate  of  Thm.  3;  in  fact  we  get  the  rate  achievable  with  1  term  (M=0)  in  the 


Fig.  1 1 .  Error  in  W  ’ 4 


-  norm  for  N 
and  2 ) ,  ^  = 

p  **  0 . 


=  0  (M  =  1 
.1414  ,  and 


€ 


Our  choice  of  tpQ,  <p^ ,  and  <P2  was  not  ideal;  it  made  some  algebra 

handier,  out  for  larger  N  one  sould  need  to  find  i|>  out  of  Thm.  1  and  2  to 

maintain  an  easily  invertible  A(w)  as  well  as  ensure  that  w.  -*•  0  fast  as 

J 

e  +  0  . 

Computations  were  also  carried  out  for  other  types  of  functions  F 
(see(1))  exhibiting  much  the  same  behavior  (although  a  less  extensive  analysis 
was  done,  see  [14])  and  also  for  the  linear  case  (F  is  the  identity)  where  it 
was  possible  to  compare  the  CPU  time  spent  with  that  of  another  solver:  FEARS 
(Finite  Element  Adaptive  Research  Solver),  cf.  [23],  solving  (4)  directly,  and 
our  method  compared  favorably  (about  30  times  faster;  it  should  be  noted, 
though,  that  FEARS  is  a  more  general  purpose  solver). 


We  conclude  that  our  approach  yields  an  efficient  numerical  procedure  for 
solving  (7)  for  c  small,  and  smooth  compatible  boundary  data  as  well  as  for 


M 

zy. 


i S 


m 


m 


i 


1 
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e  not  small  or  in  the  presence  of  some  layer  in  the  solution. 

More  precisely,  we  would  like  to  emphasize  the  following  conclusions: 

1 .  There  is  for  domains  like  ours  and  possibly  archetype  domains  a  large 
savings  in  computational  effort  compared  to  standard  nondimensional 
reduction  solvers. 

2.  The  approach  is  a  general  one,  which  no  doubt  can  be  adapted  to  other 
"stress-strain  laws"  as  expressed  through  the  function  F  in  (1)  as  well 
as  other  types  of  boundary  conditions. 

3-  The  selection  of  basis  functions  {^.}  is  robust  in  at  least  three  ways. 

<1 

a.  Even  though  the  selection  was  motivated  by  small  values  of  the  half¬ 
thickness  of  the  domain  (e)  ,  it  works  quite  well  for  e  not  small 
(e  =  1  e.g.)  without  having  to  use  an  excessive  number  (N)  of 
basis  functions. 

b.  Even  in  the  presence  of  nonsmooth  or  noncompatible  loads  giving  rise 

to  layers  in  the  solution,  the  method  works  well  by  picking  these 

layers  up  in  the  coefficient  functions  { v . }  to  the  basis  func- 

J 

tions.  These  v .  -*•  0  fast  as  j/e  *  ®  and  thus  rough  input  data 
J 

are  taken  care  of  at  fixed  e  through  increasing  N  at  least  near 
the  layer. 

c.  Through  the  parameter  p  we  have  introduced  three  asumptotic  ranges 
of  loads:  zero  (p<0),  finite  (p«0),  and  infinite  (p>0)  limit  load. 

We  showed  that  one  should  select  the  same  basis  functions  in  each  of 
the  three  cases.  Also  computationally  the  method  works  equally  well 
in  all  three  cases. 


ijh’M'vSri 


The  method  lends  itself  to  adaptivity  or  the  "feed-back"  approach  in  at 
least  two  ways. 

a.  At  fixed  N  ,  the  ODE  solver  used  is  actually  already  adaptive,  see 
Babuska  and  Majer  [21]  and  Majer  [22]. 

b.  We  would  like  to  propose  the  following  generalization  of  the  method 
in  this  paper:  let  the  choice  of  N  be  made  locally  in  each  subin¬ 
terval  of  some  given  mesh  on  the  basis  of  local  error  estimators.  We 
conjecture  that  the  computational  effort  thus  involved  would  be  far 
less  than  that  involved  in  "feed-back"  or  adaptive  methods  for 
standard  solvers  since  the  dimension  of  the  problem  has  been  reduced 
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